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1 Introduction 

Holographic superconductors are curved geometries which exhibit spontaneous breaking of 
a gauge symmetry and are dual to field theories which correspondingly exhibit spontaneous 
breaking of a global symmetry. The dual field theories are typically large ./V gauge theories 
with relativistic conformal invariance. The symmetry breaking order parameter is the expec- 
tation value of a composite operator in the field theory which is a singlet under all the gauge 
symmetries of the field theory, but transforms under the global symmetries that are gauged 
in the bulkj^] In holographic superconductors, the broken symmetry phase has a finite charge 

1 The global symmetry of the field theory which the order parameter breaks can in principle be weakly 
gauged. This is analogous to BCS theory, where phonon-electron interactions are the crucial dynamics of 
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density. This distinguishes holographic superconductors from vacuum configurations with 
spontaneous symmetry breaking. Holographic superconductors whose order parameter is a 
scalar under Lorentz transformations are called s-wave superconductors. These were first 
discussed in [TJ [2] . Holographic superconductors with a vector order parameter are called 
p-wave superconductors. Early studies include [31 HI [3]. 

Consider as an example a dual field theory which is four- dimensional and includes chiral 
fermions X a and R a , both transforming in the complexified adjoint of an SU(N) gauge group, 
and both carrying charge — e under the U(l) symmetry that gets spontaneously broken. 
Then the symmetry breaking order parameter could be (trA a A a ) in the case of an s-wave 
holographic superconductor, or (tr X a o^J^P) in the case of a p-wave superconductor, where 
a 1 is the first Pauli matrix. There is no lattice structure in the dual field theory; thus the 
breaking of rotational invariance by the p-wave condensate is spontaneous. The holographic 
treatment of both the s-wave and p-wave cases is typically at the level of mean field theory in 
that one solves classical equations of motion in the bulk without inquiring about the role of 
fluctuations. This is justified on the field theory side if one is restricting attention to leading 
order effects in a large N expansion, where N is the rank of the gauge groupj^] Notable 
reviews of holographic superconductors include EH IH] • 

One useful probe of the electronic structure of high-T c superconductors is the angle 
resolved photoemission experiments (ARPES) which essentially rely on the photoelectric ef- 
fect and measure the energy of the electrons emitted from the sample. ARPES experiments 
demonstrate some interesting properties of high T c superconductors. These properties in- 
clude a gap with d x 2_ y 2 symmetry, deformed Dirac cones whose apexes are the nodes of the 
superconducting gap, Fermi arcs whose zero temperature limits are the nodes in the gap, 
and a peak-dip-hump structure of the emission intensity as a function of frequency at fixed 
wave- number. A review of ARPES measurements can be found in |12j . 

Holographic superconductors differ from high-T c superconductors in several respects: 
Most notably, they are symmetry-breaking states of large N gauge theories, and they have 
no underlying lattice structure. Nevertheless, it is interesting to ask whether or not the 
properties of the fermionic spectral function in high T c superconductors are shared by their 
holographic counterparts. In [131 El E3 , following [T6l IT71 IT8] , an analysis of fermion cor- 
relation functions was carried out for the holographic s-wave superconductors introduced in 

symmetry breaking, and the U(l) of electromagnetism can be treated as a global symmetry for most purposes: 
For example, conductivity is computed using a two-point function of the Noether current corresponding to 
the U(l) symmetry before it is gauged. 

2 In O |8] non mean field theory behavior was exhibited by arranging for a particular bulk lagrangian 
for the scalar field. 
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[H [2] • In this work, after elucidating some details of the phase diagram of the holographic 
p-wave superconductors in AdS4, we discuss some properties of its fermionic correlation func- 
tions. Along the way we introduce a variant of the holographic p-wave construction, based 
on an 5*0(4) gauge group. 

This work is organized as follows. In section [2] we review the construction of a p-wave 
superconductor in AdS4 and study its phase diagram. An interesting property of the p-wave 
superconductor, first observed in [19] , is that the bulk geometry is a domain wall interpolating 
between infrared and ultraviolet limits which are each AdS4 with an asymptotically flat 
gauge connection^] In section [3] we discuss in detail scalar, fermion, and vector correlation 
functions in the asymptotic regions of such backgrounds, starting from a general gauge 
group. In section [4] we explain how to compute the fermion spectral function in the full, 
zero-temperature, domain-wall, p-wave background and described some of its properties. 
The numerical results of the computation can be found in section [5] In section [6] we discuss 
the results in the context of ARPES experiments on high T c superconductors. The reader 
familiar with the p-wave superconductor background, and interested only in the features of 
the spectral function should have a look at sections |3.5[ [5] and [6j The reader unfamiliar with 
computations of correlators in the context of AdS/CFT may find section [3] useful. 

2 The p- wave holographic superconductor 

The simplest example of a non-abelian holographic superconductor is the p-wave supercon- 
ductor first introduced in [4J (following earlier work (3j) and further studied in [5] [19] and 
[20] [2T| |22| [23] . The bulk action for the p-wave superconductor is given by 

S = jj'x^g (r+^-I^F 2 ^) ■ (!) 
In what follows we will set L = 1. Here F^ u is the field strength of an SU (2) gauge potential: 

F^v = dpA v - d u A fl - ig Y M [A^, A v ] , (2) 

and 

A, = AX , (3) 

3 To be precise, the gauge connection has a field strength whose stress tensor becomes insignificant near the 
boundary compared to the negative cosmological constant. It is in this sense that the ultraviolet geometry 
can be approximated by AdS4 with a flat gauge connection. 
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and r a = \o a i where a a are the Pauli matrices. Consider a configuration in which the gauge 
field takes the form: 

A = $(r)r 3 dt + W(r)r x dx. (4) 

Based on the general arguments in J3JH], one expects that apart from the Reissner-Nordstrom 
solution for which W = 0, there exist other solutions where W ^ 0, corresponding to a non- 
zero expectation value of the boundary current J*. This expectation value spontaneously 
breaks both the £77(2) gauge symmetry and rotational invariance. Symmetry breaking solu- 
tions of this type have been studied in the limit where #ym - > °o m \M an d i n a conjectured 
zero temperature configuration in [19] . In [22], five dimensional, non-zero temperature and 
finite gauge coupling geometries were studied. The purpose of this section is to fill a gap in 
the literature by studying the AdS4 p-wave superconductor geometry for finite temperature 
and coupling. Through this numerical study we will confirm that the domain wall geometries 
of [32] are indeed the zero-temperature limits of symmetry-breaking black holes. 
Parametrizing the line element by 

ds 2 = -r 2 7 (r)e-* (r) dt 2 + + r 2 {c{r) 2 dx 2 + dy 2 ) , (5) 

r A j{r) 

we find that the equations of motion for the gauge field and the metric are 

l -rW ( 2ri - - r X ' 7 + 4ry) + gXM^g + ^ w „ = Q 

2 \ c J r 2r y 



2 

1 2 _ (2d , 4\ g\ u W 2 ® 2 _„ n 
-r 2 7$' — + X ' + - - YM 2 2 + r 7$ = 

2 \ c r J r i c A 

r 2 1C " rc'(rj' + 87) g 2 M W 2 $ 2 e* , / W' 2 \ , 1 VJR/2 



+ 7 --r-=- - 3 - ri - -e x & +3 = 



2c 4r 4 c 2 7 ' V 4c 2 J '4 

We have omitted an additional equation of motion which is automatically satisfied once the 
gauge fields and metric components solve (J6]). In the limit where #ym — > 00, the matter 
content of the theory decouples from gravity and the equations of motion reduce to gauge 
fields in an AdS4-Schwarzschild or AdS4 background. This is the probe limit, initially studied 
in [I], following earlier work [2] on a similar limit of the holographic Abelian Higgs model. 
In the following subsection we will solve the equations of motion numerically and discuss 
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some of the features of the solution. We will revisit the probe limit in section 



2.2 



2.1 Numerics and phase diagram 

The equations for W, <E> and c are second order while the equations for 7 and x are nrs t 
order. Thus, to obtain a solution, we need to specify eight integration constants. In the deep 
infrared (IR), located at r = 0, we set $ — > 0, W — > finite, 7 — > finite and c — > finite. (If 

IR IR IR IR 

we are looking for finite temperature solutions then we should require that 7 vanishes at the 
horizon r = t~h, and if we are looking for zero temperature solutions we should require that 7 
is finite in the deep IR). Near the boundary, located at r — > 00, we require that W — > and 
that $ — > /i, n being the chemical potential of the boundary theory. The remaining two 
integration constants can be thought of as the values of x an d c at the boundary, and these 
can be gauged to and 1, respectively, by rescaling the t and x coordinates. In practice, 
we've looked for solutions by using a standard shooting algorithm from the horizon to the 
boundary. 

To analyze the stability of the solutions with If ^ we compute the boundary theory 
grand canonical potential Q (per unit volume) of these configurations and compare it to the 
grand canonical potential of the Reissner-Nordstrom black hole. The computation of Q is 
carried out by computing the on-shell Euclidean action Q = —TSe/V. We refer the reader 
to [22] for the details of a similar computation in AdSs. We find that solutions with a non- 
vanishing condensate, W 7^ 0, are stable only below a critical temperature T c , which varies 
with the charge #ym- Above T c the canonical potential for the AdS Reissner-Nordstrom black 
hole is lower, and it is the preferred solution (see figure [I]). We find that the phase transition 
from the RN solution to the condensed solution is second order if ^ym > 1.14 ± 0.01, and is 
first order if #ym < 1.14 ±0.01. When #ym < 0.710 ±0.001 the condensed solution no longer 
exist s|f] 

In [19] it was conjectured that the zero temperature limit of the condensed phase is a 
domain wall geometry similar to the one described in [24J: The infrared and ultraviolet 
geometries are both asymptotically AdS but with different speeds of light. The condensate 
W, the gauge field $ and the metric components x an d c interpolate between their infrared 
values Wir > 0, cjr > 1, xm an d $ir = in the infrared to their UV values Wyy = 0, 
cuv — 1; Xuv — 1 an d $uv = > 0. 7 approaches 1 both in the UV and in the IR. 
Note that since c differs in the UV and IR, the appropriate speeds of light differ in the x 



4 With 20 digits of working precision, the lowest value of q for which Mathematical NDSolve algorithm 
could obtain a solution was q — 0.7103. 
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Figure 1: The difference between the grand canonical potential for the condensed phase and 
the uncondensed phase, AQ. The red dots show the location of the critical temperature, and 
the black dots the location of the spinodal points. The left plot corresponds to gyyi = 0.79 
and the right one to #ym = 3. 
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Figure 2: The phase diagram for the p-wave superconductor. The blue curve indicates a 
second order transition and the red line a first order phase transition. The dashed black 
lines are spinodal curves and the red dot marks a tricritical point. 
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Figure 3: Plots of the metric components 7 and c at zero temperature as a function of 
the the charge. As the charge approaches its critical value, the black hole hair gets pushed 
further into the IR, conforming to an extremal Reissner-Nordstrom black hole in the UV. 
The dashed black line corresponds to the extremal Reissner-Nordstrom solution. The dashed 
vertical red line on the right plot signifies the location of the extremal RN horizon where 
7 = 0. 

and y directions. This anisotropy is a new and distinctive feature of the domain wall of [19] . 
Another interesting feature is that the AdS radius is the same in the ultraviolet and infrared: 
indeed, both limits are just empty AdS 4 with a flat SU (2) gauge connection. 

By following the branch of symmetry-breaking solutions from the region where W is 
perturbatively small to the region where it is larger, we were able to verify that the domain 
wall geometries described in the previous paragraph are indeed the zero temperature limits of 
solutions with regular, finite-temperature horizons. See figure |4j From the phase diagram of 
the SU (2) superconductor, depicted in figure [3j and by comparing the domain wall geometry 
of the condensed phase with subsequently small #ym to that of the extremal RN solution 
(figure [3]), it appears that crossing #ym = 0.710±0.001 at T = results in a phase transition 
from the condensed solution to the extremal RN solution. 

2.2 The probe limit 

As discussed earlier, in the probe limit the matter content of the theory decouples from 
the metric and the equations of motion ^ reduce to the Yang-Mills equation in a fixed 
background geometry. For convenience we rescale the gauge field to get rid of all the factors 
of <7ym m the equations of motion: that is, $ — >■ $/(?ym and W —> W/gyu- This has the same 
effect as setting #ym = 1 in the Yang-Mills equations — but we should bear in mind that the 
probe approximation is justified precisely by taking $ym large. An additional simplification 
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Figure 4: Plots of 7 as a function of the temperature for q = 0.8. The curves are color 
coded according to the temperature of the solution. The dashed line shows the domain wall 
solution and the dotted lines correspond to superheated or supercooled solutions. 



is possible at zero temperature: we recall that AdS4 is conformally flat, and that the classical 
Yang-Mills equations in four dimensions are conformally invariant. Explicitly, if the metric 
is expressed as 

L 2 

ds 2 = —{-dt 2 + dx 2 + dy 2 + dz 2 ) , (7) 
z 2 

then systematically dropping the overall prefactor L 2 / z 2 has no effect on the Yang-Mills 
equations. With the ansatz Q, they take the form 

$" = W 2 <$> W" = -<5> 2 W , (8) 

where primes denote d/dz. The boundary conditions appropriate for describing the type of 
domain wall solution we are interested in are 

$ and W -> W m at z = 

(9) 

$—>•// and W -> as z -> +00 , 

where Wx& and \x are finite. Even though one can find one conserved charge for this system 
(namely the Hamiltonian associated with radial translations), the equations of (JsJ) do 
not appear to be integrable. A closely related system studied in is known to exhibit 



S 



strongly chaotic behavior. 

We find it convenient to change variables and make the field redefinitions as follows: 

$(z) = W m $(C) W(z) = W m W(C) , (10) 

where 

C = e -^Co, (11) 

and Co is a constant yet to be determined. Large z now corresponds to ( = 0, and we can 
reformulate the boundary value problem by requiring that at C = 0, 

W(0) = 1 $(0) = $'(0) = 1 W'(0) = . (12) 

The UV is located at C = Co an d Co is the smallest C for which 

W(C) = 0. (13) 

The equation for X = W + i$ is 

C 2 d 2 ( X + Cd c X + ht{{X*) 2 -X 2 ) = (14) 
which we can solve by Taylor expanding around C = 0. Defining 

oo 

x = ]Tra n o\ (15) 



n=0 



we find that the a n 's satisfy: 



[V 1 ] k 



n 2 a r . 



^ a n -(2k+i)a2k+i- P a P (16) 

k=0 p=0 



where [(n — l)/2] means the largest integer that is smaller or equal to (n — l)/2. Equation 



(16) can be solved recursively once we are given ao = 1 and oc\ = 1. The first few terms in 



the expansion are 

w(o + mo = i + k - -c 2 - -c 3 + — c 4 + — c 5 - — c 6 - • • • (17) 

and we note that the quantities 4 n n! n!! a n are positive integers at least to n — 100. Once 
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Figure 5: The approximate functions $(C) and W(() evaluated from (15) with n = 40, with 
the positive root of W shown. 



again, by considering the first 100 terms of the series, we approximate the radius of conver- 
gence of the Taylor series expansions for X to be 3.38. Between and 3.38, W has a zero 
(meaning the phase of X is in/ 2) that appears to be unique. It is at 

Co ^2.5918. (18) 

See Fig. [5] 

One may now extract fi = $(Co)W / ir ~ 1.9285Wtr. Note that the range of ( corresponding 
to < z < oo is Co > C > 0- It is convenient that the radius of convergence of the Taylor series 



(17) is significantly larger than Co : because of this, we can get uniformly good accuracy for $ 



and W on the interval of interest from the Taylor series expansion. In practice we terminated 



the series (17) at n = 100. In section 5 we will approximate this T = background with a 



sharp domain wall solution depicted in figure [5] 



3 Two-point functions from conformally flat backgrounds 

With zero-temperature and finite-temperature backgrounds in hand, an obvious next step is 
to compute Green's functions associated with interesting physical phenomena, such as con- 
ductivity and photoemission. Such computations typically rely heavily on numerics, because 
even the coefficients in the differential equations to be solved are known only numerically. It 
turns out, however, that some of the features of the final answer can be understood qualita- 
tively in terms of the Green's functions one would get for the ultraviolet and infrared AdS4 
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geometries alone. These asymptotic forms can be obtained analytically. This section is de- 
voted to a general study of such Green's functions for more general gauge groups and even 
for slightly more general gravitational backgrounds than AdS4. The backgrounds we study 
are conformally flat, and we require that the gauge connection is also flat: that is, there is 
no Yang-Mills field strength. In the context of AdS/CFT, the flat Yang-Mills connection 
in the bulk implies that the field theory lagrangian on the boundary has been explicitly 
deformed by components of conserved currents associated with continuous symmetries, and 
that all the deformations commute both with each other and the original Hamiltonian. As 
explained earlier, our motivation for carrying out this analysis is to get a better understand- 
ing of the infrared and ultraviolet behavior of the domain wall geometries of the non-abelian 
holographic superconductors described in the previous subsection. Our analysis is, however, 
more general. We will see that we can understand not only the fermion correlators, but also 
the UV and IR limits of scalar and vector correlation functions for a generic non-abelian 
gauge group. 

The calculations in this section, though slightly abstract when expressed in terms of 
an arbitrary semi-simple Lie algebra and arbitrary unitary representations thereof, are in 
essence formal elaborations of the simplest calculations possible in AdS/CFT. We have even 
avoided the use of curved geometries by restricting attention to the conformally coupled 
scalar, the massless Dirac fermion, and non-abelian gauge fields with a Yang-Mills action. 
The equations we solve can be understood as describing the free propagation of massless fields 
in flat space — with a few mild modifications dictated by the bulk gauge symmetries. Thus, 
this section is relatively self-contained and may serve as a useful introduction to holographic 
Green's functions for formally inclined readers. 

Readers wishing to pass over the technical details but understand the qualitative features 



may wish to skip straight to section 3.5 



3.1 The action 

Consider the following action for gravity, a Yang-Mills field of a simple gauge group with Lie 
algebra g, a complex scalar field £ in the unitary representation of g, a Dirac fermion \1/ 
in the unitary representation r# of g and a real singlet scalar 0: 




(19) 
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Ordinarily there is a prefactor 1/2k 2 = 1/IQhGn multiplying 5", but we drop this factor. We 



normalize the trace tr in (19) so that tr t a t b = \5 ab ) where the t a are generators of g. Since 



the gauge group is simple, the field strength is given by 

Ffxu = d^Ay - d v A„ - ig Y m [Ap, A v ) , (20) 



and we have defined 



D^ = (d fl - i g YM A a t a r jZ 
Dfl = (d, + \u^T PJL - ig YM Alt% ) * 



(21) 



Here are the generators of g in the representation r, u^— is the spin connection, and 
underlined Greek indices are tangent space indices. When there is no risk of confusion, we 
will omit the subscript on representation matrices like : so > f° r example, t a H = i" £. We 
define gamma matrices so that T- is anti-hermitian, while the other are hermitian, and so 
that {r^, T^} = 2r]^ where jf^ = diag{-l, 1, . . . , 1}. Also, we define * = ^TK The terms 
in the action Sbd y includes boundary terms which do not affect the equations of motion but 
render the on-shell action finite. 



There are several simple ways in which we could generalize ( 19 ) — but typically will not for 
the purposes of this section. First, we could generalize to gauge group that are semi-simple or 
that contain abelian factors. The main difference is that then one has multiple independent 
gauge couplings. Second, we could introduce an explicit mass term for E, and/or other S 
dependence in V. Third, we could introduce a Dirac mass for and in the case of a real 
representation for fermions, also a Majorana mass. And fourth, depending on the choice of 
representations and r$, we could introduce some Yukawa interactions. 



The action (19) is similar to what one typically finds for the gauged supergravity trun- 
cations of the near-horizon dynamics of D- and M-brane constructions. It was chosen for 
purposes of producing simple illustrations of the Green's function computations that we 
want to do. The biggest difference between ( Jl9| and supergravity actions is that in super- 
gravity, there are spin-3/2 fields (the gravitini) in addition to g^ v , A^, S, and 0. Further 
differences are that in gauged supergravity, one would typically redefine the metric through 
a scalar-dependent Weyl transformation so that the Einstein-Hilbert term has no prefac- 
tor; and that in supergravity, the kinetic terms of the other fields would typically involve 
interesting functions of the scalars. 
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3.2 Conformally flat backgrounds 

The simplest class of solutions to (19) is 



ds 2 = e (z) 2 (-dt 2 + dx 2 + dy 2 + dz 2 ) 4> = <f)(z) , (22) 

with all other fields set to 0. For example, if V(<p) = — 6/L 2 , then 

e (z) = L/z <f)(z) = (23) 



provide a solution to the equations of motion. An solution (23) is the Poincare patch of 



AdS4, dual to a conformal field theory (CFT) in 2 + 1 dimensions. In the coordinate system 



used in (22) and (23), z runs from 0, which is the boundary of AdS 4 , to infinity. Typically 
the CFT is a large N gauge theory, and the large N approximation justifies the classical 
treatment of the bulk dynamics. The large N gauge groups of the boundary theory have 
nothing to do with the bulk gauge group g. Rather, g is the algebra of continuous global 
symmetries of the boundary theory. We will denote the Noether currents associated with 
these symmetries by J^, where a is an index for the adjoint of g and m runs over the three 
boundary directions. 



There are more general solutions to the equations of motion following from (19). In what 



follows we will be interested in solutions of the form (19) where eo and </> are defined for 
z > 0, and where eo is a monotonically decreasing function of z which diverges at z — 0. 
Such solutions can frequently be related to vacua of conformal field theories deformed by an 
operator dual to 0. When such a relation exists, the solutions are termed holographic 
renormalization group (RG) flows. True to this name, these geometries break conformal 
invariance. However, they manifestly preserve the Lorentz invariance of R 2,1 . They also 
preserve the full gauge invariance under g, because <p is (by assumption) a gauge singlet. 



What we want particularly to note is that the solution (22) (even when it's not AdS 4 ) 
is unaltered upon the introduction of non-zero gauge field A^, provided the field strength 
Ffu, vanishes. More particularly, we are interested in configurations where A z = and A m is 
constant for m = 0, 1, 2. In order to be a flat connection, A m must take values in a Cartan 
subalgebra h of g: 

A m = A a m t a eh. (24) 
Through an appropriate choice of basis, we can insist that all the t a occurring in the sum 



(24) belong to h. We will refer to the flat gauge connection (24) as a Wilson line. This 
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would be more strictly appropriate if we compactified the x and y directions into a torus 
and formed the integrals J A over closed curves 7 wrapping the torus. 



The gauge fields (24) explicitly break some or all of the gauge and Lorentz invariance, 



while preserving translation invariance. In the context of AdS/CFT, if the lagrangian of the 



field theory dual to the holographic RG flow (22) is C , then the lagrangian with the gauge 



fields (24) turned on is 



C = Co + A a m J mo 



(25) 



Two-point functions of operators which are singlets under g, such as and the stress tensor 
T mn , respect the Lorentz invariance of the background (22), at least to leading order in N, 



because the free propagation of singlet fields, such as perturbations of <fi and g^ u , do not 



respond to the Wilson lines (24). On the other hand, correlation functions of the operators 



Oy, and Oq, dual to £ and and also of the Noether currents J^, must be sensitive to the 
Wilson line. The main goal of this whole section is to describe that dependence. 

The key to making the calculation of the two-point functions of Os, Oq,, and tractable 
is that the action simplifies greatly upon the conformal transformations 



a 
Co 



3/2 " 



(26) 



The parts of the action (19) involving A„, S, and \l/ can now be expressed as 



flat 



d 4 x 



M 



-tr F 
2 tt > 



d 3 xiipr_ip 



(27) 



where we have now written 5bd y for the fermions explicitly in terms of the projection matrix 
r_, and we have defined 

1 xT- 

r ± = . (28) 



M is now the z > part of R 3,1 , and dM denotes the surface z — 0, which is just R 2,1 . 
In (28), and in the rest of this section, we use the flat metric r]^ v = diag{ — 1, 1, 1, 1} in 



all formulas: thus, for example, tr F^ u = •q^^rf 11 ' 2 tr F^ lU1 F^ 2V2 . Because there is now no 
distinction between "curved space" indices /x and tangent space indices /1, we simply use /1. 
Because the metric is flat, the spin connection vanishes. So 



D„a = {d„ - ig YU Ap a )a 
= (d„ - ig Y MAlt a ) 4 . 



(29) 
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3.3 Scalar two-point functions 

We define 

G*{x) = -i9(t)([O^t,x),OU0M) = J J^E^'* * 



(30) 



Here T denotes time-ordering, x = (t,x), k = (u,k), and k ■ x = —ut + k ■ x. We use 
(A) = tr pA where p is the density matrix of the state under consideration. When the 
temperature and chemical potentials are zero this reduces to the vacuum expectation value 



(0|A|0). The Green's functions (30) can be extracted from solutions to the wave equation 
for cr, as we now explain in some detail. 
The scalar equation of motion is 

D fl D> 1 a = 0. (31) 

Because of the translation invariance in the (t, x) directions, we can cast solutions in the 
form 

a{t,x,z) = e ik - x a(z). (32) 
A straightforward calculation shows that 

D,a = e lk - x D,a (33) 

where 

An = iK m D z = d z (34) 

K m = k m -g YM A a m t a . (35) 

{d 2 z - K m K m )a = , (36) 

so the solution is a = e ±Kz v, where K = \/K m K m and v is an arbitrary vector in r^. Recall 
that t a is in the Cartan subalgebra h of g and we are working in a basis where all the t a are 
diagonal. Thus, the K m are diagonal, and one can define non-polynomial functions of the 
K m . Indeed, let v\ be the weight space associated with weight vector A a , i.e., t a v\ = \ a v\ 



and 



(31) can be recast as 



15 



for all t a G h. Then K m v\ = (k m — k mt \)v\ where 

k m + = gYMA a m \ a . (37) 

When defining K = y/K m K m we need to specify which branch of the square root to 
pick. To start, let the spatial part k of k be large enough so that k m — k m) \ is space- 
like for all weights A. Then the action of K on vectors in the A eigenspace is simply 
Kv\ = \J{k m — k m ^\)(k m — k™)v\ = \f(k — k\) 2 v\, where we pick the plus sign on the 
square root. Thus, in this purely spacelike case, where K 2 is a positive definite matrix act- 



ing on rs, we choose K also to be positive definite. Then the solutions to (36) which decay 
far from the boundary take the form 



a = e~ Kz v . (38) 



Still assuming that K 2 is positive definite, one may extract the Green's function in the 
following manner. The on-shell action for the scalar is 

Son-shell = - / d 4 x\D^a\ 2 = - I d 4 x [d^ D"a) - (7 t D^D"a] = [ d 3 xa j d z a, (39) 

J M ' JM ' JdM 

where in the second step we used the equations of motion, and in the third step we used 
Stokes' Theorem and remembered that D z = d z . The basic premise of AdS/CFT is that the 
on-shell action is the generating functional for Green's functions of the boundary theory. In 



the current instance, this means that 



^on-sheii = W 2 [a hdy , a j hdy ] = ~ I d 3 x 1 d 3 x 2 aJ 3dy (xi)G s (xi - x 2 )(r hdy (x 2 ) , (40) 

JdM 

where Gs is a two-point function of 0%, and o"bd y is simply a evaluated at the boundary 



z — 0. The alert reader will note that we have not specified whether (44) is a retarded or 



time-ordered Green's function. At the moment we do not need to because the Hermitian 
parts of G|? and coincide, and we are doing the computation for k m such that K 2 
is positive definite. To recover the full Green's function we will use an appropriate pole 
passing prescription to go to more general K 2 . 

Now consider the following linear combination of solutions to the equation of motion for 

a: 

a = ^e ik - x e- Kz Vl + i 2 e iq ' x e~ Qz v 2 , (41) 
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where £1 and £2 are complex numbers, and Q is defined from the momentum q m as K is from 
k m . The on-shell action is now sesquilinear in £1 and £2- Plugging the boundary limit of the 



specific form (41) into the last expression in (40), one obtains immediately 



d 2 Wo 



(42) 



On the other hand, plugging (41) into the last expression in (39) leads to 



d 2 S< 



on— shell 



d 3 xvle- iq - x {-K)e tk - x v 1 



ik-x„ 



Equating (42) and (43) leads to 



Gv(k) = K 



(43) 



(44) 



This is almost our final result. The full retarded or Feynman Green's function can be 



obtained from (44) by an appropriate pole passing prescription: 



00 — > oj + ie 
u —> oj(1 + ie) 



to obtain G R 
to obtain G F . 



(45) 



The subtleties of recovering the imaginary part of Green's functions from a supergravity ac- 
tion which is real have received significant attention [271 EH EHl |3Q] . Here let us simply make 



the well-known observation that with the prescription (45) for retarded Green's functions 



the scalar wave-functions are infalling: on eigenspaces where K m acts as a timelike vector, 



K is a negative imaginary number, meaning that e 



-Kz„ 



e ipz v with p > 0. In other words, 



the number flux of scalar quanta is away from the boundary, falling into the bulk. 



Already from the simple expression (44) we can reach one of the main conclusions of this 



section: The spectral measure is supported in the timelike regions of light-cones in fc-space, 
one through k m \ for each weight A in the representation r^. To see this explicitly, we can 
simply note that Gs(k) acts as 



G s , x {k) = ^{k 



kx) 2 



(46) 



on the weight space associated with A. This part of the Green's function has an imaginary 
part precisely when k — k\ is timelike, i.e. in the light-cone passing through k\. The spectral 
measure is essentially this imaginary part. 

We have not concerned ourselves with the overall normalization of the action. If we had, 
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we would have found an overall prefactor on which scales as N 3 / 2 in constructions based 
on iV coincident M2-branes. 

Giving the scalar E an explicit mass term m, instead of the conformal coupling — ||S| 2 
implies that the dual operator O^, has dimension A where A(A — 3) = m 2 L 2 and L is the 
radius of curvature of AdS4. The calculation of the two-point functions proceeds almost 
identically to what we laid out above, except that it has to be done in AdS4 rather than 
flat space, because the mass term is not conformally invariant. The radial parts of the 
wave-functions (without any conformal rescaling) come out proportional to z 3 ^ 2 'K^_ 3 / 2 (Kz) 



instead of ze~ Kz as in (38). The pole passing prescription (45) will convert the modified 
Bessel function Ka-3/2 into an appropriate Hankel function when K m acts as a timelike 
vector. When A ^ n + 3/2 with n an integer, the result for the two-point function is 



Gs{k) = c&K 2A ~ 3 , where ca is a fc-independent factor. Evidently, (44) corresponds to the 
case A = 2. When A = n+ 3/2 the Green's function takes the form C/\K 2A ~ 3 log K. The 
limit 2A — 3 — > n coincides with this form of the Green's function since ca diverges in this 
limit. 

3.4 Spinor and vector two-point functions 

Our results for the scalar Green's functions can be summarized in the following simple 
terms. With no gauge field present in the bulk, conformal invariance dictates GY,(k) = 
ca& 2A ~ 3 - For A = 2, the same expression can be deduced more generally for conformally 
flat backgrounds. The effects of the gauge field are simply to replace k with K, where 



K m = k m - gyuA^t 11 as in (|35j). 

In light of these results, it is reasonable to guess that Green's functions for spinor op- 
erators and conserved currents in the presence of a flat connection follow from similar re- 
placements k m — > K m . More explicitly, the Green's functions we have in mind are defined 
as follows: 

f rl 3 k 

G*(x) = -ie(t)({O 9 {t,2),Ol{0,0)}) = / 7^e*"G«(*) 

J ( 27r ) (47) 
r d 3 k 1 ' 

G%?(x) = -ie(t)([J^(t,x),J b n (0,0)}) = / iA-J k - x G^{k). 



, (2tt) 3 

Time-ordered Green's functions would be defined similarly. We usually omit the adjoint 
indices a, b. may have any dimension A, but because is conserved, its dimension 
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must be 2. Our expectation is that 

>y m K ( K 2 n — K K \ ab 

where /a arid Soo are constants that we do not propose to track at this stage, and the ie 



prescription (45) is implied. Our notation, essentially following (30], is to represent the 



four-dimensional gamma matrices as 





r 2 = „ J , (49) 

yy-" u y 
where 

7* = ^2 7 1 = °"i 7 2 = °"3 
and a a denotes the Paul matrices. We decompose the Dirac spinor ip as 

*=(*+]■ («) 




The operator transforms as a two-component spinor, like ?/>_. We sometimes consider 
the Dirac conjugate of a two component spinor: ip± = ^±7*- The Green's function Gy(k) 
has the same spinor structure and group representation content as the bilinear if)_ip_. 

Just as in the case of scalar Green's functions, if U- is a two-component spinor which 
belongs to a weight space of with weight vector A then the action of on m_ is of the 
form Gq,(k)v,- = Gq, t \(k)u_, where 

G/I \ J- 7 ("771 k m \) x /rr>\ 

*,a(*) = ~/a (fc _ fcA)4 _ 2 k 7 f , (52) 



and fc^ is defined as in (37), only using weights of r# not r s . Likewise, if e n is a polarization 
vector which also belongs to a root space of g with root vector a, then G mn {k)e n = G mn ^ a {k)e n 
where 

. — y // \ (^ ^a) Vmn y^m k ma ^(k n k na ) , , 

"mn.a^J = s oo ; , == j (53) 

- k a y 

and fc TOjQi = ,gYM^4m aa - Evidently, the spectral measures of and G mn are supported inside 
light cones, just as for G^, but each light cone has its apex at the the momentum k\ or k a 
associated with a weight vector A or a root vector a, as appropriate. 



Let's sketch a derivation of the expression in (48) for Gy in the case A = 3/2, correspond- 
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ing to a massless fermion, again focusing first on the situation where K 2 is positive definite. 
This is the case where we can do all the calculations in flat space, where the equation of 
motion is 

™^ = 0. (54) 



This equation implies D^D^ip = 0, which is identical to (31). So the allowed solutions are 
of the same form: 

if; = e ik - x e~ Kz u . (55) 



Here u is a four-component spinor transforming in the representation r#. Plugging (55) into 



(54) provides a constraint on u: 



{iT m K m - T Z K) u = . 



(56) 



Decomposing u as in (49), one can rewrite (56) as 



.l m K m 



(57) 



To compute the fermionic Green's function we plug a linear combination of solutions to 
the Dirac equation, 

^ = ^ e ik - x e- Kz Ul + £ 2 e iq - x e- Qz u 2 , (58) 



into the on-shell action, which is simply the boundary term in (27): 



Son-shell = - / d 3 xiipT_ip = -i£*£i(27r) 3 <5 3 (g - fc)w 2 +«i_ + • • • ■ 



(59) 



In (|59|) we omitted terms involving different combinations of the & and their conjugates. 

(60) 



Using (57) we find 



d 2 S, 



on— shell 



On the other hand, essentially by definition 



2vr) V(g - k)u % 



l m K n 



K 



-u 1+ 



d 2 W 2 



-(2n) 3 S 3 (q - k^+G^k^u^ 



(61) 



Comparing (60) and (61) leads immediately to 



= — j^r 



(62) 
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Continuing to momenta where K 2 < can be done using (45). 

For the gauge field, the story is essentially the same. We start by perturbing the gauge 
field: 

A = A^ t + a. (63) 
The field strength is F = ^(d — igA) 2 = f + 0(a 2 ), where 

ffu, = D^a v - D v a^ , (64) 

where by convention D = d — igAft at . The linearized equation of motion is 

^ = 0. (65) 

It is most straightforward to proceed in a gauge where a z = 0. Then one can show directly 



from (65), assuming as usual that K 2 is positive definite, that a m = e tk ^ xm e Kz e m if K m e r 



0, and that a m = e lkmXm e m if Kn j^" - e m = e n . Here e n is a constant polarization vector. The 
transverse modes, where K m e m = 0, are physical and the others are pure gauge. If we define 

K K m 

P n m = S n m - , (66) 

then we may compactly write the allowed solutions as 

a = e ik x e- KPz e , (67) 



where Lorentz indices for the boundary directions are now implied. Plugging (67) into the 
quadratic on-shell action 



Son-shcii = J d 3 x tr a m f zm , (68) 
one finds more or less immediately that 

G mn (k) = KP mn , (69) 



which indeed has the form indicated in (48). 



We will not delve deeply into conductivity calculations in this paper, but it is worth 



noting that (69) contains a geometrical explanation of the hard-gap phenomenon remarked 
on in [19] . To understand this connection, first recall that the gauge potential in the infrared 
copy of AdS4 is (in the notation of [19]) A = Bqt 1 dx. So the Cartan subalgebra h is the 
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one generated by r 1 . The conductivity studied in [T5] is the one related to a perturbation 
5A = e~ luJt a(z) t 3 dy. Now, G^ n (k) acts in the adjoint representation, so it can be denoted 
Gmn b (.k). The conductivity of interest is 

a yy ^) = -v a v b G^ y a \uj,0,0), (70) 

where v a is the unit vector corresponding to the generator r 3 . Now we need to find a 



convenient basis in which to work out the right-hand-side of (70). The obvious basis for 
the adjoint representation, for our purposes, is the one in which r 1 acts diagonally. The 
eigenvalues of its adjoint action are a = ±1 and 0: these are the roots a of the adjoint 
representation of SU(2), and they correspond to the c-number Green's functions G^ ina 
defined in ( [53] ) . In the basis for the adjoint representation where (1,0,0) corresponds to 
a = 1, (0, 1,0) corresponds to a = —1, and (0,0, 1) corresponds to a = 0, the unit vector 
corresponding to the generator r 3 is v = (1, 1, 0)/v2. Plugging this expression for v into 



(70), we have 

i 



a yy>T s(u;) = ^ [G£ |+ (w,0,0)+G£ > 4w,0,0)] . (71) 



Now we can explain the geometric origin of the gap. From the explicit expression (53) we 
see that G yy + is real except inside the light cone whose apex is at k m ^ + = gyM-Bo^m- The 
distance along the u axis one must go in order to cross into this light cone is gyM-Bo- Likewise, 
Gy y _ is real except inside a lightcone which one crosses into along the u axis at the same 
value of lo. G yyQ has no gap, but it is not involved in a yy ^{uj). So the gap is A = g YM B . 

What we are neglecting is that there is asymmetric warping of time and space in the 
domain wall geometry of [19] . The t coordinate picks up a warp factor e~ X0//2 in the infrared 
AdS4 geometry relative to the ultraviolet, while the x coordinate picks up a factor cq. Corre- 
spondingly, the wave- vector k m+ at the apex of the light cone gets scaled by l/c because it 
points in the x direction, and frequency gets scaled by e X0//2 . Altogether, introducing these 
rescalings into the expression A = gyM-Bo for the gap leads to 

A = ^^ymBo (?2) 
c 

A should match the gap found in [19], and it does. (Recall that their q is our ^ym-) What 
we learn in addition is that if conductivity could be measured at finite wave-number, the 
gap would decrease as one goes toward the apex of one of the lightcones, and can be reduced 
to zero if the wave-number is chosen to lie at the apex. Similar analysis seems possible for 
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<7 XXT 3 , but we have not fully considered the consequences of mixing with the graviton and 
with timelike components of the gauge potential in this case. 

It may seem striking that we are able to evaluate the gap by considering only the infrared 
limit of the geometry. This actually makes sense because there is a continuum contribution 
to the spectral measure of G^ b only when the corresponding gauge-boson wave-function has 
some infalling component in the infrared geometry; otherwise the Green's function is purely 
real, with no dissipative part, except for ^-function localized contributions to the spectral 
measure like the one that signals infinite DC conductivity. The same reasoning, essentially, 
is behind the evaluation of the gap in [TH]: the Schrodinger potential plateaus in the infrared 
at an energy that determines the gap. 



3.5 Summary and an example 

The upshot of the calculations in this section is that the Green's functions of operators dual 
to the conformally coupled scalar S, the massless Dirac fermion and the non-abelian 



gauge fields in (19), are simple functions of modified momenta 

K m = k m -g YM A a m t a . (73) 

In particular, they all involve fractional powers of K 2 = i] mn K m K n . K 2 is always hermitian, 
but positive definite only for k outside all the shifted light -cones {k m k m ^\){k m — 0. 

As explained earlier, k m ^\ is given by 

k m ,x = 9YMA a m X a , (74) 

where A is a weight vector for the representation in which the operator whose Green's func- 
tion we are considering transforms. The spectral measure of these Green's functions comes 
entirely from the region where K 2 fails to be positive definite, i.e. the region where k — k\ 
is timelike for some weight A. 

It is interesting that the Green's functions have a pure power law form not only for an 
AdS4 bulk, but for any conformally flat bulk geometry which is asymptotically AdS4. This is 
a consequence of the invariance of the relevant part of the classical supergravity action under 
conformal transformations. Thus, power-law correlators, like those one sees in conformal 
theories, arise at leading order in a large N limit for special operators in backgrounds that 
are dual to non-conformal theories. At subleading orders in N, loop corrections in the 
bulk enter in, and one cannot expect the exact power-law behavior of the boundary theory 
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correlators to persist unless there is conformal symmetry. All the backgrounds treatable by 
the methods of this section have exact boost invariance, 5*0(2,1), as well as translational 
invariance. These symmetries would be preserved at all loop orders in the bulk, meaning all 
orders in N in the boundary theory. 

As an example, consider g = so(4) = su(2)a © su(2)b, where the generators for the two 
su(2) factors are labeled r\ and r%. The standard Cartan subalgebra is generated by t\ 
and r%. Let the fermions \l/ transform in the 4 = 2a x 2b of so(4). Consider now the flat 
connection 

A = W (dxT A + dyrl) , (75) 
where W is a constant. Assuming that for X = A, B the eigenvalues of each r£ acting on the 



doublet of su(2)x are ±1/2, the images (74) of the weight vectors of the vector representation 
of so(4) are 



k 



gyMW 



SlS 2 



Sih + s 2 k 2 ) 



(76) 



where s« = ±1 and k\ are unit vectors in momentum space. Thus the apexes of the Dirac 
cones pass through points at the corner of a square whose sides are aligned with the k x and 



The spectral measure of the retarded version of the Green's function (62) is 



p(co,k) = -rmtrG|(&) = 2\cj\ ^ 



9{cu 2 



\k k SlS2 \ 2 ) 



Si=±l 



\k — h I 2 



(77) 



where we have used (62) and (76). The most conspicuous feature of (77) is that p(u,k) is 
supported in the union of four Dirac cones, whose apexes are the vectors k SlS2 . Plots of the 
Dirac cone and of p(u>, k) as a function of cu for fixed k are shown in figure ki| By adding a 



mass to the fermion, we can adjust the power in the denominator of (77) from a square root 



to 2 — A, where A is the dimension of the dual fermionic operator in the field theory: see 
PI). 



4 Properties of the spectral function 

We are interested in studying the Fourier space retarded Green's function G R (k m ) for 
fermions in the doublet of SU(2) for the zero temperature p-wave superconductors dis- 



cussed in section |2| Our conventions for G can be found in (47). Note that G has both 
spinor and gauge indices and a total of 16 a priori independent components. Here we will 
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Figure 6: Features of the spectral function for the operator dual to a massless Dirac fermion 
in the 4 of so(4), in a conformally flat background with the flat so(4) connection (75). This 
is distinct from the spectral function in the p-wave superconductor to be discussed in the 
next section. Top: The Dirac cones inside which the spectral measure is supported. Bottom: 
Two representative plots of p(co, k) as a function of u for fixed values of k. 
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be interested in the spectral measure, p, denned by p = — ImtrG , where the trace is over 
both the gauge indices and the spinor indices. This quantity is gauge invariant, and it is 
interesting to consider because it gives a measure of the number of eigenstates of the theory 
that couple to the fermionic operator, with long lived excitations manifesting themselves as 
(^-functions peaks. In what follows we explain in some detail how to compute the spectral 



function (section 4.1 ), how it relates to normal modes of bulk fermions (section 4.2) and how 



positivity of p is guaranteed from a bulk point of view (section |4.3 ). The reader interested 
in the final numerical results for p may skip directly to section [5j 

4.1 Green's functions from a gravity dual 

In section [3] we computed the retarded Green's function for a configuration with a flat 
connection. In this section we consider the fermion Green's function for the zero temperature 
p-wave superconductor discussed in section [2] where the connection is flat only in the IR and 
UV. The action for the background is given by ([I]) and we use the notation in (|4]|5]) to describe 
the metric and SU(2) gauge fields. We introduce the spin-1/2 field \I/ which transforms in 
the doublet of 577(2), whose action is given by 

S {eimion = -i [ d 4 x^^D^-i [ d 3 x^/^r^T^ , (78) 

JM JdM 



where is given in (21). Our conventions for the T matrices are largely as in section 3.4 
but because we now work in curved spacetime, we must be careful to distinguish between 
curved and flat indices. Because the matrices 7" 1 are needed to describe the physics of the 
boundary theory, which is defined on a flat background, we persist in defining them as in 



(|50|). In place of (|49|) we employ 

■ ' * = -r 4 = r I . (79) 




The second term in (78) is a boundary term that does not affect the equations of motion 
but gives the only nonzero contribution to the on-shell action. We treat \l/ in the probe 
approximation and do not allow it to back-react on the geometry. The equation of motion 
for \l/ is D^ty = 0. To write it explicitly, it is convenient to exploit translation invariance in 
the x m directions and write \l/ in the form 

\S?(x m , r ) = (-gg rr y\ e - iult+ik * x+ik y y ilj(r) . (80) 
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Note that (80) reduces to (26) for a conformally flat metric. As in section [3| it is useful to 



split ip into two chiral spinors 



4 




The equations of motion for %p± take the form: 



X 

e 2 



(-00 - 9ym$ r 3 ) 7* + - - g YM W r 1 ) 7* + k y ~f y 



^7> 2 



X 

e 2 



,_ (-w - 3ym$ t 3 ) 7* + - - gymWr 1 ) 7* + ^7* 
-V7 c 



"0H 








(82) 



where the prime denotes a derivative with respect to r and 7 M are 2 + 1 dimensional bound- 



ary theory 7 matrices related to the bulk T matrices through (49). Note that ip± have 



both a boundary spinor index and an 577(2) doublet index. We suppressed both types of 
indices in (82) but it should be clear how they are contracted. For instance (rfr 1 ^-^)^^ = 
ryt w) r i(«*)^,u) ( 6 ) ) where (i), (J), . . . = 1, 2 are spinor indices and (a), (6), . . . = 1, 2 are S77(2) 
doublet indices. 



The asymptotic solution to (82) in the infrared is 



Ui - V ie - K+/r + V 2 e~ K - /r + V 3 e K+/r + V A e K -' 

2 



2c 



IR 



Here c — >■ cir and similarly W — > Win an d X — * Xm- The spinor integration constants V\ 

IR IR, IR 

and V2 satisfy t 1 V,i = ±^Vi with the minus sign for % = 1,3 and the plus sign for i = 2,4. 



They must also satisfy a linear constraint that follows from (56). The equations k± = 



define the two IR Dirac cones that we expect for SU(2) doublet fermions. As discussed 
in section [3] if we are outside one of the Dirac cones, i.e, if for instance n 2 + > 0, then the 



corresponding term in (83) should be regular provided we pick the positive sign for the 



square root in the definition of k + . If we are inside one of the Dirac cones, the correct 
boundary conditions for computing the retarded Green's function are obtained by making 
the substitution uj — > u + ie and demanding regularity. In practice, this prescription amounts 
to keeping uj real and picking the sign of the square root such that the imaginary part has 
the same sign as uj. This is equivalent to requiring that the solutions are infalling in the IR. 
Thus, we set 

V 3 = Vi = . (84) 
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The asymptotic behavior of the fermions in the UV takes the form 



^ UV = Qie- X+/r + Q 2 e x+/r + Q 3 e 

1 \ 2 P _\ f 

A ± = | (a;±-p YM $uvJ ( . 2 

where c — )• cuv an d a similar definition for %uv an d $uv- The are constant spinors 
such that r 3 Qi = ±^Qi with the plus sign for i — 3,4 and the minus for i — 1,2. Like Vi 




and V2, each Qt satisfies a linear constraint that follows from (56). As opposed to the IR 



asymptotics (83), in the UV there is no restriction on the choice of the sign of the square 
root in the definition of \±. The equations \± = define two UV Dirac cones, different from 
the IR Dirac cones k± = 0. The IR Dirac cones are shifted in the k x direction because in 
the IR only A x is nonzero. The UV Dirac cones are shifted in the uj direction since in the 
UV only A t is nonzero. Furthermore, since c and \ go to different constants in the UV and 
IR, the UV and IR Dirac cones are stretched by different factors in the u and k x directions: 
that is, they are characterized by different, anisotropic speeds of light. 

The Green's function for the operator dual to ^ can be found from the UV behavior of 
the fermions following the prescription given in [311 E2] and more recently in [30] . For the 
purpose of computing the Green's function we rewrite the UV asymptotics in the form 

^ = A + B- + O fl 
^ = D + C- + 2 

where A, B, C and D are arbitrary chiral spinors with an SU(2) doublet index which are 



1 > 86) 



linearly related to the Qi. The linear constraints that follow from (56) can be used to solve 



for B and C in terms of D and A. Since the equations of motion are linear, requiring that 



the solutions are regular in the IR, (84), imposes a linear relation between A and D, 



As explained in section [3J to obtain G R we take the functional derivative of the on-shell 
action. The result is 

G R = -iM-i 1 , (88) 
where we suppressed the spinor and SU(2) indices. The derivation is very similar to the one 
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for uncharged fermions, and we refer the reader to p3U] for details. 



4.2 Normal modes 

We are particularly interested in the regions in which p ^ 0. Consider the IR asymptotics for 



the fermions, (83). If n\ > or k?_ > (i.e., we are inside the IR light-cone) then ip + can be 
chosen to be everywhere real and consequently, p = except perhaps for poles of Retr(Gpt) 
which, together with the ie prescription corresponds to p having 5-function support. In what 



follows, we argue that a divergence of Gr corresponds to a solution to (82) which is a normal 
mode. 



We define a normal mode to be a solution of (82) which is regular in the IR and its near 



boundary expansion takes the form (86) with A = 0. This is a natural definition of a normal 
mode since A is the leading term in the UV asymptotic expansion of if). We can not require 
D to vanish as well since we are not free to choose D once we fix A. Normal modes will only 
occur for special values of k m and Vi (i.e., the values of k m for which there are normal modes 
will form a codimension one surface). The coefficients A and D in the UV expansion of ijj 
and V\ and Vi in the IR expansion of ip are both integration constants of the same solution 
to a linear differential equation. Hence, 

A = UV D = VV, (89) 

where V should be understood as a four component object that contains the four independent 



integration constants in the V,. Using (87) the four by four matrices U and V are related to 
M. through 

M = VU-\ (90) 

where some care must be taken in interpreting the index structure. Since V, A and D are 
finite the entries of U and V can not diverge. From here it follows that if G R has diverging 
entries then det hi = 0; the only way that Ai (and hence G R ) can diverge is if U is not 
invertible, i.e., the solution is a normal mode. 

By the same reasons as argued in [15] , we expect normal modes only outside the IR Dirac 
cones. Let us briefly repeat this reasoning. If k\ > then, as discussed earlier, ip + can be 
chosen to be everywhere real and the condition det U = is a real equation for which we 
can expect to find solutions for appropriate u, k x and k y . On the other hand, if k\ < 0, 
ip + cannot be chosen real and det U = will be a complex equation whose solutions are 
expected to involve complex uj. Such a solution would correspond to a quasinormal mode 
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rather than a normal mode. 

There is also a somewhat weaker argument for expecting the normal modes to be inside 



at least one of the UV Dirac cones, (A+ < or A < 0). According to (85), if A+ > and 
A 2 > then the generic UV solution will grow exponentially and matching such a solution 
to the IR asymptotics seems improbable. If at least one of A+ and A 2 is negative, then 
there is at least a two dimensional subspace of oscillatory solutions and there is no such 
obstruction. Our numerics indicate that the outermost UV light cone is indeed where the 
surface of normal modes ends. We conclude that the normal modes exist in a "preferred 
region" 

k\ > and k 2 _ > and (a+ < or A 2 < o) . (91) 

This preferred region is bounded provided the UV Dirac cones have a narrower opening angle 
in both the k x and k y directions than the IR Dirac cones. This is certainly true of the SU(2) 
backgrounds we have constructed. It may be possible to demonstrate in general that the 



region (91) is compact starting from an appropriate positive energy condition. 



4.3 Positivity of the spectral measure 



Unitarity requires that the spectral measure p = — Im tr Gr is nonnegative for all real k m in 
any field theory, as can be shown using a spectral decomposition. Instead of explaining the 
spectral decomposition argument (which is standard) we will show in this section how the 
positivity of the spectral measure follows from a computation in the gravity dual. 
Consider the current 

= -fp 1 * . (92) 
which is conserved in the sense that V M J M = provided the equations of motion that follow 



from (78) are obeyed. Not surprisingly, this conserved current is associated to the gauge 



invariance of the action (78). Here V is a covariant derivative. From (80), we can write 



V M J" = as Q'{r) 



0, where 



Q(r) = y/^gJ r = -ipr^j = + = ~ 2Re ■ 

where in the last line we used = — Vv0-- Using (86) we see that 



(93) 



Q = -2 Re (A a V L> a ) . 



(94) 
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Combining (94) with (88), it follows that 

g = 2Im(4 )7 *G Jl W (6)7 *AW) . (95) 

Since, A is an arbitrary spinor and 7* is invertible, Q < implies that (G — is a 
nonpositive definite matrix and therefore, p = — ImtrG is everywhere nonnegative. 



It remains to show that with our IR asymptotics (83), Q < 0. To this end we work out 



the constraint on the Vi that follows from (57). If we split the V{ into chiral spinors, denoted 
by V^, the constraint can be written as 

Vr = i°^V+ V 2 - = i—V+ (96) 

where 

a ± = -e 2 ^7 H 7- + ^7^. (97) 

2cir 

There are now several different case we must consider. Suppose first that we are outside 
both IR Dirac cones. If we take the to be real and choose the Majorana basis for the 7 



matrices (making them real matrices) it follows from (96) that the V[~ are pure imaginary. 
Therefore, ip+ip- is also pure imaginary and Q is identically zero leading us to conclude that 
p = (up to, perhaps, a codimension one surface of normal modes.) Next, suppose that 
k 2 + < and take Vj> = for simplicity. If we choose real then V{~ will also be real and Q 
is no longer zero. It is given by 



Q = j^-j- Re I y+Va+VT 1 • (98) 

It is now a simple exercise to compute the eigenvalues of the matrix 7*a + and show that it 
is positive (negative) definite for u > (u < 0). Therefore, Q is always negative for k\ < 
and V2 = 0. We can show in the same way that Q is negative for h?_ < and V\ = 0. Since 
Q is given by a quadratic form, it follows that Q > 0. 



5 Evaluating the spectral function 



5.1 Numerical results 



To obtain the spectral function p = — Im tr G R we solved (|82j) numerically with the initial 



conditions given in (83) for some small but nonzero 7*j. We then extracted A and D from 
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Figure 7: The normal mode surface for ^ym = 3. The IR and UV Dirac cones are always 
drawn red and purple respectively. The normal mode surface is contained in the preferred 
region (91) and is shown as black dots. 



the behavior of the numerical solution at large r. This procedure was repeated for the four 



linearly independent choices of the Vi in (83). We then used (87) to obtain M. and (88) to 
obtain Gr. To find the normal modes we looked for zeros of det U. 

Let us now discuss our numerical results. The spectral measure p(k m ) is symmetric 
separately under the sign flips uj —> -co, k x — > —k x , and k y —> —k y . Since we are restricting 
our attention to massless fermions the only free parameter is #ym, the Yang-Mills coupling 
constant. Most of the features discussed in this section seem to be generic and independent 
of <7ym> with the obvious restriction that (Jym is at least large enough so that the T = 
background we are considering exists. It seems possible that if we take gyM large enough 
there could be several disconnected surfaces of normal mode, as was found in [15] for the 
s-wave background, but we were unable to check whether this happens as our numerics are 
not reliable for all of the "preferred region" for large gyu- F° r 9ym as large as 10, there is 
only one distinct surface. 

The surface of normal modes we obtained numerically can be roughly described as a 
truncated deformed version of the IR Dirac cones, as can be seen in Fig. [7j This surface 
does not extend indefinitely and is contained in the preferred region (91). It is perhaps more 
interesting to consider the shape of the normal mode surface near u = 0. The reason for 
this is that, in our conventions where $(0) = 0, a fermion with oj = has energy equal to 
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the chemical potential^] For the p-wave superconductor, Fig. [7] clearly shows as uj — > the 
surface of normal modes approaches the IR Dirac cones and that therefore the Fermi surface 
is given by two isolated points: the apexes of the IR Dirac cones. 

To understand the behavior of p inside the IR Dirac cones it is useful to plot p(k m ) as 
a function of uj for fixed k x , k y . Some of these plots are shown in Fig. [8j Outside the IR 
Dirac cones (the boundaries of which are shown as dashed red lines), p vanishes except for, 
perhaps, a 5-function peak that indicates the presence of a normal mode. Inside the IR Dirac 
cones, p is positive and smooth except for the intersection of two Dirac cones where there is 
a kink. 

For k m close to the apex of a Dirac cone, it is natural to expect that we will recover the 



behavior extracted from the infrared asymptotics. This behavior was discussion in section 3.5 



for the SO (4) case, and it is essentially the same in the 577(2) case except in that there are 
only two Dirac cones rather than four. It is clear from plot (3) in Fig. [8] that p goes to a 
constant as uj — > when k is at the apex of the Dirac cone. Plot (1) in Fig. [8] reveals that 
there is almost a power-law singularity as soon as one crosses into the Dirac cone close to its 
apex. One does not see this dramatic near-singularity when crossing into a Dirac cone far 
from its apex. 

To gain a better understanding of the spectral measure we provide some plots of constant 
uj slices in figure |9j At small uj we find that the normal modes surround each of the IR light 
cones. For the particular case of qym = 3 depicted in the top left corner of figure [9] the 
normal modes terminate on the light-cone leaving a "scar" in its interior. This is also seen 



more clearly in figure 10 This ridge inside the IR light-cone suggests that the normal mode 
has turned into a relatively long-lived quasinormal mode, or resonance. 

Going back to figure [9j we observe that as uj is increased the shape of the surface of normal 
modes becomes more asymmetric until, eventually, the surface of normal modes around each 
of the light cones intersect (figure [9] (3)). At large uj the normal modes arrange themselves 
into an inner an outer closed surface. This inner surface disappears once the IR Dirac cones 
start overlapping. The outer surface will also disappear, eventually, and this happens as 
soon as the outermost UV Dirac cone is inside the IR Dirac cones. 



5 If we wanted to treat VP beyond the probe approximation, the natural construction is to fill all the 
excitations with energy below the chemical potential, i.e, all the normal modes with oj < 0. Once this is 
done, the Fermi surface will be given by the intersection of the normal mode surface with the ui = plane. 
But we don't understand how to systematically treat the back-reaction of the fcrmions on the geometry and 
the gauge field. 
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Figure 8: The preferred region and the spectral measure p as a function of a; for fixed k x , k y . 
The IR Dirac cones are shown in red, and UV Dirac cones in purple in the top figure. The 
blue lines are constant k x and k y lines and are the axes of the plots on the bottom. The 
limits of the IR and UV Dirac cones are marked in these plots as vertical dashed lines that 
are red and purple respectively. All the plots are for g YM = 8. 



34 



tuff* 




(1) (2) 




-10 1 -10 1 

kjjfi kjfi 

Figure 9: The spectral measure p for u/p = 0.45, 0.52, 0.53, 0.58 and #ym = 3. The UV 
and IR Dirac cones are shown in the top figure as red and purple surfaces respectively. The 
blue planes are the constant u planes for which p is plotted in the bottom figure. In these 
plots, outside the the IR Dirac cones the surface of normal modes is shown as a black curve. 
This is overlaid with a density plot of the spectral measure p where black corresponds to 
large values of p. 
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Figure 10: A contour plot of the spectral function p on a constant u/p = 0.93 slice for 
9ym — 8- The red and purple curves mark the location of the UV and IR cones respectively. 
The black line outside the IR light cone shows the location of the normal modes. The spectral 
function p develops a ridge between the points where the normal modes terminate on the 
light-cone. 



5.2 Fermion correlators in the sudden approximation 

The numerical work reported so far in this section turned up some interesting features in 
fermionic Green's functions: multiple thresholds associated with overlapping Dirac cones, 
normal modes in a preferred region, long-lived quasi-normal modes, and some signs of recov- 
ering the infrared approximations near the apex of a Dirac cone. We would like, if possible, 
some analytic approximation that exposes these features more clearly, so that we are not so 
dependent upon numerical integration of classical equations of motion. The obvious place 
to start is the probe approximation discussed in section 2.2| This approximation relies on 
having g YM big so that the gauge field doesn't back-react appreciably on the AdS 4 geome- 



try. Recall that in section |2.2| we rescaled fields as needed to eliminate explicit dependence 
of the equations of motion on #ym, and that at T = we were able to replace the AdS 4 
geometry by the z > half of R 3,1 . For the sake of simplicity we will continue to restrict 
to a massless Dirac fermion, which (after a suitable rescaling) we denote if). It is dual to a 
spinorial operator Oq,- We do not restrict if) at this stage to be a doublet of SU(2): it will 
become apparent that most of our discussion can be carried through straightforwardly for 
a domain wall configuration based on any semi-simple gauge group, with if) transforming in 
any representation of it. 

The two-point functions of are controlled by solutions to the linear equation of motion 
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for i/j: 



r»{d li -iA li )ii> = o. 



(99) 



Because we are working in the z > half of R 3,1 , we do not need to distinguish between 
and T^. is a domain wall solution to the flat-space Yang-Mills equations, so for the 
SU (2) case it is takes the form 



SU(2): A m dx m = §(z)dtT 6 + W(z)dx 1 r 1 



(100) 



where <3> and W satisfy the equations fl8J). For 50(4), an interesting solution generalizing 



(75) is 



SO(4) : A m dx m = $(z) dt (t% + t%) + W(z) (dx 1 r\ + dx 2 t% 



(101) 



with the same functions $(z) and W(z) as in the SO (4) case. Backgrounds based on more 
complicated gauge groups could also be constructed. All we require for the discussion of the 
next few paragraphs is that A m depends only on z, and A z = (a gauge choice). 



Solutions to (99) can be cast in the form 



ip(x, z) = e ■ x ip(z) 



(102) 



where, as usual, k m = (—lj, k) and x m = (t,x). If we define 



(103) 



then (99) becomes 



(d z + tT z T m K m )^ = 0, 
whose solutions can be formally expressed as 



ip(z) = P |exp 



-1 I dz'T z T m K m {z') 




1>(0) 



(104) 



(105) 



where P denotes ordering non-commuting matrices so that those coming from larger values 
of z go to the left. The allowed solutions satisfy boundary conditions in the infrared (large 
z) which can be compactly expressed as 



ip oc e 



-K m z 



U . 



(106) 
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Figure 11: The functions $ and W, with Wxr = 1, together with the sudden approximations 
^sudden and Madden discussed in the main text. 



where restrictions on u come only from solving (104), and we define 



= lim K m (z) 

z— >co 



K™ = lim K m (z) 

z— >oo 



IR m 



n UV ll m ) 



(107) 



and the ie prescription (45) is implied. 



It is hard to go further than (105) without some additional approximation because the 



path-ordered exponential is hard to compute. The domain wall structure suggests an ob- 



vious approximation, illustrated in figure 11 Let's replace A m — >■ A;^ dden , where A s ^ dden is 



piecewise constant and piecewise flat, going straight from the UV flat connection to the IR 
flat connection at a special value z* of z such that 



dzA 



sudden 



dz A+ . 



'108) 



A straightforward calculation based on the function $(C) introduced in section 
that in the SU(2) and 50(4) cases described by ( flOOj ) and ( floTj ), 



2.2 



shows 



W m z, « 1.2058. 



(109) 



From now on we will set Wm = 1 for simplicity. It is easy to solve (99) with A m replaced 
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by A s £ dden : the solution is 

V'sudden(^) — 



zT z T m K}™ 



^sudden(O) for < Z < Z* 

^sudden(O) for Z* < Z . 



(110) 



Now let's see how to implement the boundary condition (106). Because -^sudden satisfies the 



Dirac equation (with the replacement A m — > A™ dden ), we have 



d z *l> S udden(z) = -irT m i^^ sudden (z) for all z > z*. 



(HI) 



On the other hand, d z Sudden (z) = — ifrR,VWdeii(<z) for large z because of (106). Since the 
iC? commute with one another, it must be that 



3*$Buddffli(s) = -^lRVWden(z) for all z > z* 



;ii2) 



Comparing (111) and (112), we arrive at 



^sudden(O) = 



where 



P = (K m - iT z T m K l *)e 



lK\-iz,T z T n K^ v 



It will be useful to express 



P = q + T z T m q m + T m T n q n 



(113) 

(114) 
(115) 



where 



q = K m cosh(z*if 



UV, 



K™cot&(z.K uv ) + K m K i 



r 



uv 



uvj 



uv 



(116) 



n - ^IR^UV Sinll (^ KuV ) 



K 



uv 



Note that q, q m , and g mn have no spinor structure: they are purely gauge-theoretic quantities. 



In (115), the first term on the right hand side is implicitly multiplied by the unit matrix in 



spinor space. 



The equation (113) determines the fermionic Green's function in the sudden approxima- 



39 



tion, as we now explain. If we use our usual basis for gamma matrices, (49), and express 

"^sudden (0) = 




(117) 



then (113) can be recast as 




where 



r 




H 


\p + 


p-j 





q _|_ 7 m 7 n g mn 



-i m q m 



J m q m 



q + l m l n q m n , 



:ii8i 



(119) 



Comparing to (86) and using (87), we arrive at 



GWdenW = iP+-P++l l = tPZ-P-+l l 



(120) 



A more explicit form, based on the middle expression in (120), is 

G sudden (k) = -«(7 m <?m)" 1 (g + 7 m 7 n, ?mn)7* 



:i2r 



We further define p su dden(&) = — I m ^r G^, dden (/c). We note that trG^ udden (fc) and hence 
Psudden(^) can in principle be found in closed form as functions of k m , /i, Wjb,, and z*. In 
practice, the closed-form expressions for q, q m , and q mn in terms of k m , /i, W^, and 
are already quite complicated, and we were unable to find a closed-form expression for the 
inverse matrix (7 m (?m) _1 that was sufficiently compact to be useful. 

Although we did not succeed in finding a simple enough closed form expression for the 



spectral measure to record here, the expression (121) is simple enough to expose most of the 
qualitative features of the analytic structure of G su dden(^)- First we claim that branch cuts 
in q and q m , as functions of k m , arise only when K\r has a branch cut, while q mn has no 
branch cuts at all. To demonstrate this claim, we observe that cosh^-fTuv) arid smh (^"uv) 
are analytic functions of K^y, which in turn is a quadratic expression in the momenta k, m \ so 



sinh(z*_R'uv) 
^UV 

and also have no branch cuts. Our first claim now follows by inspection of the 



cosh^-Kuv) arid smn ^*-^uvj have no branch cuts at all as functions of the k m . More trivially 



formulas (116) for q, q m , and q mn . 

It follows from our first claim that the continuum part of p S udden(&) is supported precisely 
where K m has a branch cut, which is to say inside the Dirac cones. This feature of the 
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spectral function has been discussed extensively in section [4] and is essentially as expected 
on intuitive grounds: The retarded Green's function can have a dissipative part iff some 
component of the fermion wave-function is infalling in the infrared, rather than exponentially 
decaying there. 

Our second claim is that, at least for ip in a real representations of g and for generic 
values oik = (k x , k y ), ^-function singularities in p su dden(w, k x , k y ) as a function of uj can only 
occur outside the IR Dirac cones. Recall that a 5-function in p su dden(w, k x , k y ) is associated 
with a pole in Retr G^ dden (u;, k x , k y ) for real uj. Thus our claim is that any such pole must 
arise outside the Dirac cones. To demonstrate the claim, first note that q, q m , and q mn never 
diverge. So the only way to get a pole is if 7 m g m is non-invertible, which is to say det i7 m g m 
vanishes. The 7 m may be chosen in a Majorana basis, where all entries are real: indeed, 
(50) is such a basis. In a real representation of the gauge group, all the and are 
real symmetric matrices, as are cosh(z*Kuv) and sinh (^*^uv) ^ ^ r j g a ^ go rea j an( ^ S y mme t r j c 

provided we are outside the Dirac cones. Thus the iq m are real matrices, and we see that 
det ij m q m must indeed be real. More precisely, for real k x and k y , det i^ m q m is a real function 
of the variable uj outside the Dirac cones. Inside the Dirac cones, Km. has an imaginary part, 
and det i7 m g m is not a real function of uj. All that we need to note now in order to complete 
our argument is that real analytic functions of a single variable generically can have zeros 
on the real axis, but general complex analytic functions do not. It is tempting to speculate 
that this argument could be extended to fermions in complex representations. Certainly 
there is intuitive reason to think that when a delta-function contribution to the spectral 
measure crosses into a continuum, it will spread out into a finite-width resonance — as we 



saw numerically in Fig. 10 



We caution the reader that the sudden approximation is not controlled in the sense of 
becoming good when some parameter is taken large or small. Usually, sudden approximations 
are justified when wave- functions are slowly varying as compared to the features of the 
underlying background that one is approximating. Optimistically, one might expect our 
sudden approximation to be good near the infrared light-cone, because then the fermion 
wave-functions are slowly varying in the region z > z*. But these wave-functions are not 



necessarily slowly varying for z < z*. Thus we regard ( 121 ) as useful in the sense of providing 
an in-principle closed-form expression that captures some of the relevant physics: namely a 
continuous part of the spectral measure inside the Dirac cones, with the possibility (at least 
on genericity grounds) of normal modes only outside the cones. 
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6 Discussion 



The starting point of our analysis is the classical action 

S = J d A x^g \ R + - i tr F* v - i^VD^j + boundary terms , (122) 

which is essentially the lagrangian of QCD coupled to gravity with a negative cosmological 
constant, except that we choose the gauge group to be SU(2) or 50(4), while the fermion 
transforms either as the doublet of SU(2) or the fundamental 4 of 5*0(4). The lagrangian 



(122) describes the bulk dynamics dual to a field theory in 2 + 1 dimensions whose continuous 



symmetries form the same group as the gauge group in (122). We treat the dynamics of 



(122) classically, which is understood to be dual to a large N approximation in the field 
theory. 

One output of our analysis is the phase diagram of superconducting black holes based 
on the SU(2) gauge group. This phase diagram is shown in Fig. |2} We demonstrated, 
largely through a numerical study, that black holes charged under the r 3 generator of SU(2) 
spontaneously break that symmetry through a p-wave condensate similar to the one originally 
studied in A conspicuous feature of the phase diagram is a tricritical point separating 
second order and first order behavior at the symmetry breaking phase transition. We also 
showed that at low temperatures, the symmetry-breaking solutions approach the AdS4-to- 
AdS4 domain wall geometries of |19j . similar to domain walls found in the Abelian Higgs 
model in [21] except for anisotropic alteration of the coordinate speed of light in the infrared. 
Our analysis is not complete in that we did not systematically study the stability of the 
symmetry-breaking solutions, and it is possible that there are other symmetry-breaking 
configurations that we missed. Thus we cannot rule out the existence of a more complicated 
phase diagram than we plotted in Fig. [5J with (for example) symmetry-breaking phases 
present below #ym = 0.710. 

A second output of our analysis is two-point functions of the fermionic operators dual 
to These two-point functions show some intriguing parallels with ARPES data on high- 
temperature superconductors. Relationships between holographic fermionic correlators and 
ARPES data were emphasized early in [T7| in the context of the non-superconducting phase, 
following earlier work [161 EH HE] • Studies in the superconducting phase include [HJEEH]; see 
also [15J. These works all focused on rotationally symmetric backgrounds. By contrast, our 
Fermi surface at T = consists of isolated points: two in the SU (2) example, and four in the 
50(4) example. Above each isolated point, a Dirac cone rises, as can be seen in Fig. [8] for 
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the SU(2) case and in Fig. [6]for the SO (4) case. In the SO (4) case, the nodes are positioned 
at 45° degrees relative to the axes along which the gauge potentials are aligned, reminding 
us of the positioning of nodes in the gap in d x 2_ y 2 superconductors. 

The structure of normal modes is also favorable to the comparison with ARPES data. 
As shown in Fig. [7| there is a normal mode slightly outside each Dirac cone. And as shown 
in Fig. [8] (plot (2) especially), the spectral measure significantly away from the tip of the 
Dirac cone exhibits a peak-dip-hump structure. The peak comes from the normal mode, 
which shows up in the spectral measure at T = as a 5-function, like an infinitely sharp 
quasi-particle. The hump comes from the continuum part of the spectral measure, which is 
entirely inside the Dirac lightcone. At the tip of the light cone (plot (2) of Fig. [8]) or just 
slightly away from it (plot (1)), there is less structure: the peak goes away or merges into 
the hump. Again we see a point of comparison with ARPES: the classic peak-dip-hump 
structure arises away from the node in the gap. 

There is a twist in our discussion of Dirac cones, normal modes, and continuum structures 
relative to the usual story based on quasi-particles, where continuum structures arise because 
of two- or three-particle states, where each particle by itself is on-shell when its momentum 
lies on the Dirac cone. In our case, the Dirac cone characterizes the edge of the continuum 
rather than the dispersion relation for the quasi-particle. We see already from formulas like 



(62) that if only the infrared dynamics are accounted for, a continuum supported inside each 
Dirac cone is the only feature of the spectral measure. There are no quasi-particles in sight in 
this infrared limit. The quasi-particles (or, at least, the normal modes of the bulk fermions) 
come from the more intricate domain wall structure of the full bulk geometry, as discussed 



in section 4J2 The dispersion relation of these quasi-particles is not perfectly linear. This is 
striking because, even in the full domain wall geometry, the continuum part of the spectral 
measure is supported in Dirac cones which are perfectly linear. In fact, close inspection 
of Fig. [7] shows that the normal modes cross into the Dirac cones in certain regions. This 
behavior is brought out better in Figs. [9] and 10 Moreover, the normal modes disappear 
altogether once one passes outside the preferred region described in (91 ). So for large enough 
bj and k there are no normal modes, and therefore no well-defined quasi-particle ^-functions 
in the spectral measure. But one still finds that the edge of the continuous part of the 
spectral measure defines a perfect Dirac cone. In short, the continuum part of the spectral 
measure is the fundamental feature, while the quasi-particle ^-function appears only under 
the right circumstances. 



We were able to produce an explicit formula ( 121 ) which captures the main features of the 
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fermionic Green's function, including the perfect Dirac cones enclosing the continuous part of 
the spectral measure and the normal mode outside the Dirac cones. Previous works, notably 
[33] . have provided analytic approximations to interesting fermionic correlators; moreover, 
the ones in [33j are based on a controlled approximation, whereas ours is not. The analytic 
forms found in [33J for the non-superconducting state rely upon the existence of an AdS2 
near-horizon region. This feature of the geometry is double-edged: While it does provide 
tractable asymptotics, it also forces the existence of non-zero entropy density at zero tem- 
perature, which remains unexplained and seems to us peculiar in a theory whose underlying 
formulation is a continuum field theory rather than a lattice. Limited analytic information 
about fermionic two-point functions is available in the backgrounds studied in [31], where 
also entropy vanishes linearly with temperature. A more powerful understanding might be 
forthcoming if one better exploited the AdS% near-horizon region of the ten-dimensional 
embedding of these backgrounds. 

There are some good reasons to be suspicious of the relevance of our setup to real high-T c 
materials with a c?-wave gap: 

• As already remarked in the introduction, the field theory dual to the AdS± bulk is a 
large N field theory formulated in the continuum rather than on the lattice. 

• The condensate in the SO (A) case is not described in terms of a spin-2 bulk field, as one 
might expect, but rather in terms of gauge potentials involving off-diagonal generators 
of 50(4). It's not at all clear that the classic phase-sensitive features of the d-wave 
gap would show up in our system. 

• Off-diagonal gauge potentials are dual to persistent currents of global symmetries in the 
boundary theory. This seems rather different from the usual language for discussing 
(i-wave superconductivity. It would be interesting to see if the notion of persistent 
currents as an order parameter might be translated into lattice language, and how it 
might interact with constraints such as those discussed in |35j . 



In addition to the 5*0(4) gauge symmetry, the field theory dual to (122) has 5*0(2, 1) 
Lorentz symmetry and also relativistic conformal invariance. These symmetries are 
largely broken by the condensates, but the fermion operator transforms as a dou- 
blet under 50(2, 1). This has no immediate analog in spin systems relevant to high-T c 
materials, where the spins are doublets under the SU (2) spin and have no further struc- 
ture under the Lorentz group in 2 + 1 dimensions. 
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The £0(4) symmetry of the Hubbard model is composed of SU (2) spin and SU(2) pseudospin . 
The definition of the latter seems to require the lattice. Doping amounts to adding 



a chemical potential for the r 3 component of pseudospin, whereas in (101) we added 
equal chemical potentials for the r 3 components of both SU (2)'s in SO (4). 

Nevertheless, the resemblance of our results for the spectral measure of fermionic Green's 
functions to the spectral properties revealed in real materials by ARPES are striking enough 
that we should inquire what underlying physics is driving it. At one level there is no puzzle: 
as soon as we note that the infrared Green's functions depend on Lorentz-invariant combi- 
nations of K m = k m — gYMA m , where k m = (-co, k) and A m is a flat connection in the gauge 



group of (122), we see that the displaced Dirac cones are just a consequence of the Lorentz 
invariance plus the eigenvalues of the non-zero components of A m . At another level, it may 
seem strangely suggestive that 5*0(4) is the symmetry group of the Hubbard model on a 
bipartite lattice, and the fermion creation and annihilation operators transform as the 4 of 
SO (4) F] Precisely this choice of gauge group and fermion representation gave us the Dirac 
cone structure reminiscent of d x 2_ y 2 pairing. Did we get approximately right answers for 
ARPES-like spectra because we have captured some correct features of the Hubbard model? 
A positive answer to this question would be fairly exciting. 



Acknowledgments 

We thank D. Haldane, D. Huse, C. Mathy, and S. Sondhi for useful discussions. This work 
was supported in part by the Department of Energy under Grant No. DE-FG02-91ER40671 
and by the NSF under award number PHY-0652782. FDR was also supported in part by 
the FCT grant SFRH/BD/30374/2006. 



References 

[1] S. S. Gubser, Breaking an Abelian gauge symmetry near a black hole horizon, Phys. 
Rev. D78 (2008) 065034, |arXiv: 0801 .2977] . 

[2] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, Building a Holographic 
Superconductor, Phys. Rev. Lett. 101 (2008) 031601, |arXiv:0803.3295| . 

6 A theory of phase competition between anti-ferromagnetic order and superconductivity has been ad- 
vanced based on an approximate SO(5) symmetry |36j ; for a review see |37j . 



45 



[3] S. S. Gubser, Colorful horizons with charge in anti-de Sitter space, Phys. Rev. Lett. 
101 (2008) 191601, |arX iv:0803.3483] . 

[4] S. S. Gubser and S. S. Pufu, The gravity dual of a p-wave superconductor, JHEP 11 
(2008) 033, larXiv:0805.2960| . 

[5] M. M. Roberts and S. A. Hartnoll, Pseudogap and time reversal breaking in a 
holographic superconductor, JHEP 08 (2008) 035, |arXiv : 080573 898 . 



[6] S. Franco, A. Garcia-Garcia, and D. Rodriguez-Gomez, A general class of holographic 



superconductors, arXiv:0906. 1214, 



[7] F. Aprile and J. G. Russo, Models of Holographic superconductivity, arXiv: 0912. 0480 



S. Franco, A. M. Garcia-Garcia, and D. Rodriguez-Gomez, A holographic approach to 



phase transitions, Phys. Rev. D81 (2010) 041901, [ar Xiv : 0911 . 1 354] . 

[9] S. A. Hartnoll, Lectures on holographic methods for condensed matter physics, Class. 
Quant. Grav. 26 (2009) 224002, [arXiv: 090373246 . 



[10] C. P. Herzog, Lectures on Holographic Superfluidity and Superconductivity, J. Phys. 



A42 (2009) 343001, |arXiv: 0904. 1975 



[11] G. T. Horowitz, Introduction to Holographic Superconductors, arXiv: 1002 . 1722 



[12] A. Damascelli, Z. Hussain, and Z.-X. Shen, Angle-resolved photoemission studies of the 
cuprate superconductors, Rev. Mod. Phys. 75 (Apr, 2003) 473-541. 

[13] J.-W. Chen, Y.-J. Kao, and W.-Y. Wen, Peak-Dip-Hump from Holographic 



Superconductivity, arXiv : 0911 . 2821 



[14] T. Faulkner, G. T. Horowitz, J. McGreevy, M. M. Roberts, and D. Vegh, 



Photoemission 'experiments' on holographic superconductors, arXiv: 0911 .3402 



[15] S. S. Gubser, F. D. Rocha, and P. Talavera, Normalizable fermion modes in a 



holographic superconductor, arXiv: 0911 .3632 



[16] S.-S. Lee, A Non-Fermi Liquid from a Charged Black Hole: A Critical Fermi Ball, 
Phys. Rev. D79 (2009) 086006, |arXiv: 0809. 3402] . 

[17] M. Cubrovic, J. Zaanen, and K. Schalm, String Theory, Quantum Phase Transitions 
and the Emergent Fermi-Liquid, Science 325 (2009) 439-444, |arXiv: 0904. 1993] . 



46 



[18] H. Liu, J. McGreevy, and D. Vegh, Non-Fermi liquids from holography, 
larXiv: 0903 .24771 



[19] P. Basu, J. He, A. Mukherjee, and H.-H. Shieh, Hard-gapped Holographic 



Superconductors, arXiv: 0911 . 4999 



[20] P. Basu, J. He, A. Mukherjee, and H.-H. Shieh, Superconductivity from D3/D7: 
Holographic Pion Superfluid, JHEP 11 (2009) 070, [arXiv: 0810 .3970] . 

[21] C. P. Herzog and S. S. Pufu, The Second Sound of SU(2), JHEP 04 (2009) 126, 
|arXiv: 0902. 0409] . 

[22] M. Amnion, J. Erdmenger, V. Grass, P. Kerner, and A. O'Bannon, On Holographic 



p-wave Superfluids with Back-reaction, arXiv: 0912. 3515 



H.-B. Zeng, Z.-Y. Fan, and H.-S. Zong, Superconducting Coherence Length and 
Magnetic Penetration Depth of a p-wave Holographic Superconductor, 
larXiv: 0912 .49281 

[24] S. S. Gubser and F. D. Rocha, The gravity dual to a quantum critical point with 
spontaneous symmetry breaking, Phys. Rev. Lett. 102 (2009) 061601, 
|arXiv: 0807. 1737| . 



[25] A. Yarom, Fourth sound of holographic superfluids, JHEP 07 (2009) 070, 



|arXiv: 0903. 1353]. 



[26] G. K. Sawidy, Yang-mills classical mechanics as a kolmogorov k-system, Phys. Lett. 
B130 (1983) 303. 

[27] D. T. Son and A. O. Starinets, Minkowski- space correlators in AdS/CFT 

correspondence: Recipe and applications, JHEP 09 (2002) 042, |hep-th/ 020505 1 1 . 

[28] C. P. Herzog and D. T. Son, Schwinger-Keldysh propagators from AdS/CFT 
correspondence, JHEP 03 (2003) 046, |hep-th/0212072| . 

[29] S. S. Gubser, S. S. Pufu, and F. D. Rocha, Bulk viscosity of strongly coupled plasmas 
with holographic duals, JHEP 08 (2008) 085, | arXiv: 0806. 0407) . 

[30] N. Iqbal and H. Liu, Real-time response in AdS/CFT with application to spinors, 
Fortsch. Phys. 57 (2009) 367-384, jarXiv: 0903 .2596]. 



47 



[31] M. Henningson and K. Sfetsos, Spinors and the AdS/CFT correspondence, Phys. Lett. 
B431 (1998) 63-68, |hep-th/980325l] . 

[32] W. Mueck and K. S. Viswanathan, Conformal field theory correlators from classical 
field theory on anti-de Sitter space. II: Vector and spinor fields, Phys. Rev. D58 



(1998) 106006, hep-th/9805145]. 



[33] T. Faulkner, H. Liu, J. McGreevy, and D. Vegh, Emergent quantum criticality, Fermi 
surfaces, and AdS2, |arXiv:0907.2694[ 



[34] S. S. Gubser and F. D. Rocha, Peculiar properties of a charged dilatonic black hole in 



AdS5, arXiv:0911.2898 



[35] F. Bloch, Flux quantization and dimensionality, Physical Review 166 (1968), no. 2 
415-423. 

[36] S.-C. Zhang, A unified theory based on so(5) symmetry of superconductivity and 
antiferromagnetism, Science 275 (1997), no. 5303 1089-1096. 

[37] E. Demler, W. Hanke, and S.-C. Zhang, so(5) theory of antiferromagnetism and 
superconductivity, Rev. Mod. Phys. 76 (Nov, 2004) 909-974. 



48 



